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A  NEW  CLASS  OF  CROSS  SECTIONS  FOR  USE  IN  ATOMIC  SCATTERING 

CAK2JIATICNS 


INTRODUCTION 

In  the  process  of  calculating  the  range,  energy  deposition  and  damage 
characteristics  of  heavy  ions  in  bulk  materials,  one  needs  sane  representation 
of  the  interactions  between  the  inpinging  ions  and  the  target  constituents. 
Because  the  target  atoms  are  displaced  and  themselves  traverse  the  target,  one 
needs  the  interactions  between  the  various  constituents.  There  is  available  a 
wide  range  of  such  potentials,  from  the  universal  ones  of  the  Thcmas-Fermi 
type  to  ones  specifically  tailored  to  reproduce  the  crystalline  properties  of 
particular  species. 

These  interactions  are  used  in  three  general  types  of  radiation  damage 
calculations.  The  most  exacting  are  the  molecular  dynamics  codes,1'2  which 
simulate  a  small  portion  of  a  crystal  in  the  acnputer.  The  response  of  the 
crystal  to  an  initial  displacement,  caused  by  an  incident  neutron  or  heavy 
ion,  is  followed  by  simultaneously  solving  the  equations  of  motion  of  all  of 
the  atoms  in  the  sanple  crystal.  More  flexible,  if  less  thorough,  are  such 
simulation  codes  as  those  of  Robinson  and  Torrens3  and  Beeler  and  Besco.4  In 
these  the  crystalline  structure  of  the  material  is  present  in  seme  fashion, 
but  the  moving  ions  are  followed  one  at  a  time.  The  molecular  dynamics  codes 
require  potentials  that  reproduce  the  major  properties  of  the  crystal.  The 
second  type  of  simulation  code  is  somewhat  less  stringent,  in  part  because  the 
potentials  are  not  allowed  to  act  for  inpact  parameters  greater  than  50-80%  of 
the  nearest  neighbor  distance. 

Finally,  there  are  the  transport  theory  calculations,  5,6,7,8,9,10  which 
assume  an  amorphous  material — no  crystal  structure  is  present.  These 
calculations  do  not  use  the  potentials  directly,  but  use  the  associated 
differential  cross  sections.  For  the  most  part,  they  make  use  of  the  of  the 
Lindhard,  Nielsen  and  Scharff  (INS)5  universal  cross  section  based  on  the 
Thanas-Fermi  interaction.  There  are  a  number  of  problems  associated  with  this 
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interaction.  It  is  known11  be  too  strong  for  large  separations.  While  its 
universal  character  makes  it  convenient  to  use,  it  does  not  accurately  reflect 
the  properties  of  more  realistic  potentials,  especially  at  large  separations . 
Hiencmenological  potentials  are  usually  made  to  vanish  for  sufficiently  large 
distances  (2nd  or  3rd  nearest  neighbor  distances) ,  whereas  the  Thanas-Fermi 
potential  is  of  infinite  range.  The  corresponding  Thanas-Fermi  cross  section 
blows  far  small  energy  transfers  and  the  total  cross  section  is  infinite, 
unless  a  cutoff  in  energy  transfer  is  imposed. 

<Xnr  goal  in  this  work  is  to  create  new  differential  cross  sections  (more 
specifically,  families  of  them)  that  both  more  accurately  reflect  the 
properties  of  realistic  potentials  and  are  as  easy  to  use  as  is  the  INS  cross 
section..  In  the  next  two  sections  we  describe  the  Thanas-Fermi  potential  and 
its  relation  to  various  phenomenological  interactions  and  outline  the  method 
that  INS  use  to  obtain  their  cross  section.  We  then  present  a  new  class  of 
cross  sections  that  reproduce  the  INS  form  for  large  energy  transfers 
(corresponding  to  small  separations  of  the  interacting  ions) ,  but  which  have 
flexible  characteristics  for  small  energy  transfers. 

Within  the  spirit  of  the  approximations  used  by  INS,  there  is  a  method, 
outlined  in  Appendix  B,  that  allows  the  potential  corresponding  to  a  given 
differential  cross  section  to  be  easily  calculated.  We  present  a  series  of 
figures  of  the  potentials  that  are  generated  from  our  new  families  of  cross 
sections.  By  comparing  a  given  phenomenological  potential  with  the  families 
thus  generated,  one  can  automatically  obtain  a  differential  cross  section 
closely  corresponding  to  that  phenomenological  potential.  Finally,  we  give  a 
number  of  examples  and  discuss  the  possible  use  of  these  new  cross  sections. 

INTERA3XMIC  POTENTIALS 

Before  considering  atonic  scattering  cross  sections,  let  us  look  at  same 
of  the  interatomic  potentials  cn  which  the  cross  sections  are  based.  Suppose 
we  have  an  atom  with  atomic  number  and  atonic  mass  A^  interacting  with 
another  atom  with  atonic  number  Z2  and  mass  A2.  Firsov12  represented  the 
interaction  between  two  such  atoms,  separated  by  a  distance  r,  as  a  screened 
Coulomb  potential  of  the  form 


The  quantity  Cp(r)  is  the  Thanas-Fermi  screening  factor,  which  is  available  in 
tabular  form.13  Die  screening  radius  a  is  given  by 


‘-[iii1]173  [7^-2]  Z'V3  =  0.8S534  ao  S'1/’ 

-  ,  Z  V2  +  Z  V*  ,* 

zl-(z12',3  +  z22/3)V2 
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where  in  expression  (3)  we  present  both  the  Firsov  version  for  Z  and  that  due 
to  INS.  This  Firsov  or  Thanas-Fermi  (as  it  is  more  cccmonly  called)  potential 
is  universal,  in  that  it  can  represent  the  interactions  between  any  pair  of 
atans,  given  only  their  atonic  numbers  and  the  separation  between  them. 

Because  of  its  simplicity,  this  potential  is  widely  used  in  atonic  scattering 
problems. 

While  the  Thanas-Fermi  potential  is  cn  reasonably  secure  theoretical 
ground  when  the  interatomic  separations  are  small,  it  is  inadequate  for  large 
separations.  The  Thanas-Fermi  screening  factor  Cp(x)  falls  off  as  x-3  for 
large  x.  Cn  both  experimental  and  theoretical  grounds  an  exponential 
falloff  is  more  realistic.  For  this  reason,  another  universal  potential,  the 
Moliere,14' 15  is  also  widely  used.  The  Moliere  potential  also  has  the 
form  (1) ,  and  its  screening  factor  reproduces  the  Thanas-Fermi  screening 
factor  out  to  reduced  separations  of  x^6.  After  that,  it  falls  off 
exponentially.  The  Firsov,  INS,  or  sane  other  suitable  screening  radius  may 
be  used  with  the  Moliere  potential. 

The  best  interaction  to  use,  of  course,  is  one  that  is  tailored, 
phenomenologically,  to  the  two  interacting  atans  in  the  solid.  There  are  many 
ways  of  obtaining  such  potentials,  but  a  ccrmonly  used  one  is  sane  version  of 
the  following  (for  the  moment,  we  assume  Z1=Z2) :  For  small  separations,  the 
Thanas-Fermi  interaction  is  used.  For  large  separations  a  form  is  used  whose 
parameters  are  adjusted  to  reproduce  seme  of  the  properties  of  a  crystal  made 
up  of  these  atans.  Finally,  these  two  portions  are  joined  smoothly  in  the 
intermediate  region.  In  Fig.  1  we  shew  such  an  interaction,  due  to  Erginsoy, 


INTERATOMIC  SEPARATION  [r  (A)] 


Fig.  1  —  Iron- iron  interatomic  potentials.  Shewn  are  the  Thomas- 
Fermi  and  the  Moliere  using  the  Firsov  screening  length  (0.0996  A) , 
the  Moliere  with  the  Torrens-Robinson  screening  length  (0.0781  A) , 
and  the  phenomenological  potential  of  Erginsoy,  et  al. 
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Vineyard  and  Englert, 2  for  alpha  (toe)  iron.  We  note  that  this  potential 
passes  through  zero  at  a  separation  somewhat  less  than  the  nearest  neighbor 
distance  (2.48  A) .  This  behavior  is  typical;  it  is  the  attractive  portion  of 
the  potential  in  the  nearest  neighbor  region  that  binds  the  crystal.  More 
sophisticated  potentials  will  oscillate  about  zero  out  to  2nd  and  3rd  nearest 
neighbor  distances;  the  magnitude  of  these  oscillations  is  of  the  order  of 


1  eV. 


Far  onrpa risen,  we  also  show  the  Thanas-Fermi  and  Moliere  potentials  for 
iron,  using  the  Firsov  value  (0.0966  A)  of  the  screening  radius.  We  see  that 
the  agreement  is  not  bad  for  small  separations,  as  we  would  expect,  but 
becomes  increasingly  worse  as  the  separation  increases,  especially  for  the 
Thanas-Fermi  interaction.  This  again  is  typical;  the  Thanas-Fermi  is  too 
large  for  large  separations.  Finally,  we  show  the  Moliere  potential  with  the 
Torrens  and  Robinson16  value  (0.0781  A)  of  the  screening  length.  Here  the 
agreement  is  quite  good  for  large  separations.  Of  course  varying  the 
screening  length  may  distort  the  agreement  at  small  separations. 

Far  the  purposes  of  transport  theory  calculations,  the  details  of  the 
potential  in  the  region  near  and  past  the  point  at  which  it  first  passes 
through  zero  are  of  no  consequence.  No  significant  amount  of  energy  is 
transferred  to  lattice  atoms  when  the  impact  parameter  lies  in  or  beyond  this 
region.  Further,  the  differential  cross  section  for  an  infinite  range 
potential  blows  up  as  the  energy  transfer  vanishes,  which  inposes  an 
unnecessary  complication  an  integral  equation  calculations.  Rather  than  using 
a  cross  section  based  on  the  infinite  range,  Thanas-Fermi  potential,  these 
calculations  would  be  both  sirrpler  and  more  realistic  if  we  started  with  a 
finite  range  potential. 


REVIEW  OF  THE  INS  METHOD 


In  computer  simulation  calculations  one  uses  the  relevant  potential 
directly  to  represent  the  interaction  between  atoms.  In  integral  equation 
calculations,  one  used  the  differential  cross  section  derived  from  that 
potential.  Suppose  that  one  particle  with  energy  E  and  impact  parameter  p  is 
moving  toward  another  (stationary)  particle.  The  center  of  mass  energy  is 


Er  ~  a2  E  /  (Ai  +  A2) 


mm 
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TOW? 


The  scattering  angle  is  given  byx/ 


«  ^  r  o  Z.  Z-  e  _  i 

9  »  r  -  2p  *  [r2[l  -  ^  C(r/a)  ]  -  p2 


where  rm-jn  is  the  largest  root  of  the  radical  in  the  integrand.  The  energy 
transferred  in  the  collision  is 


T  =  Tm  sin2  e/2 


Tm  -  r  E  =  [4  A1  A2  (A1-f-A2)“2]  E  , 


where  Tm  is  the  maximum  kinematically  allowed  energy  transfer.  Corresponding 
to  this  maximum  energy  transfer  is  the  min-imim  energy  that  can  be  carried  away 
by  the  incident  particle,  which  is  given  by 


Enin  -  E//3  -  E  (Ar-A2)2/(A±+A2)2  . 


Much  of  the  work  of  Lindhard  and  his  coworkers  is  couched  in  terms  of  the 
dimensionless  variables 


e  -  E  /  EL  »  E  [ZX  Z2  e2  (Aj+A2)  /  (a  A2)  ]-* 


t  *  e2  T/Tm  =  e2  sin2  0/2  . 


In  terms  of  these  definitions  the  differential  cross  section  is  given  by 


do  =  2  w  p  dp  ; 


INS  use  the  notation 


d a  =  ~(ir  a2  /2  )  f3/2  f(t!s)  dt  , 


so  that 


Ji 


f(t*)  =  (-4/a2)  t3/2  p(t)  dp(t)/dt  .  (13 

For  lack  of  any  standard  name,  we  will  call  f  (<p)  the  kernel  of  the  cross 
section. 

Out  of  a  desire  to  create  a  simple,  universal  cross  section,  LNS  make  two 
approximations18  to  obtain  f (t^)  for  the  Thomas-Fermi  potential:  they  replace 
Eq.  (5)  by  the  momentum  approximation19  to  the  scattering  angle,  so  that 


“  dr  [r2^2]-1/2  £ 


r  Jp 


and  they  use  the  expression 


t  *  1/4  e2  Q2 


to  replace  ©  by  t,  instead  of  the  exact  expression  (io) .  By  combining  these 
approximations  we  find  that  we  can  write 


tV2 _ ?  * 

2a 


dx  [  x2  -  (p/a)2  ]~1/2  ^  [  x"1  C(x)  ] 

.  P/a 


Eq.  (16)  provides  a  functioned,  relationship  between  the  impact  parameter  and 
the  reduced  energy  transfer  that  can  be  used  to  solve  Eqs.  (11-12)  for  f  (t'2) . 
The  advantage  of  the  INS  approximations  is  that  the  kernel  of  the  differential 
cross  section  depends  only  on  one  variable;  Eq.  (16)  is  a  relationship  between 
t  and  p/a,  not  t,  p/a  and  e. 

Given  these  considerations,  the  total  (macroscopic)  cross  section  can  be 
written 


Nd(E)  -  Nna 


dtp  <p~2  f  (0) 
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where  N  is  the  number  density  of  lattice  atans  and 


0  ~  t3*  . 


Following  IMS,  we  can  write 


2  r  Er 

N  S(E)  -  Nka  — M  d 0  f(0) 
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far  the  stopping  power  for  elastic  collisions  and 


N  W(E)  =■  Nma 


V  E.1 2  fe 
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02  f(0) 


far  the  square  fluctuation  in  energy  loss. 

The  INS  method  we  just  outlined  could  be  used  for  other  potentials,  as 
well  as  the  Thcmas-Femi.  We  will  use  a  subscript  L  to  refer  specifically  to 
the  INS  results  for  the  Thcmas-Femi  interaction.  The  quantities  fL,  SL  and 
sr&  provided  by  INS  in  tabular  fora.  A  convenient  fit  to  fL,  in  the  form 

fw(0)=  k  01/3  [1  +  (2k)2/3  08/9  ]~3/2  ,  k=l. 309  (21 

has  been  created  by  Winterbon,  Sigmund  and  Sanders  (WSS)  .8  with  this  form, 
relatively  convenient  expressions  can  be  found20,21  for  ctw,  Sw  and  Vfy;  in  some 
cases  they  are  easier  to  use  than  the  tabular  forms  of  INS. 

NEW  CROSS  SECTICNS 

In  our  discussions  of  interatomic  potentials,  we  implied  that  it  would  be 
useful  to  have  a  potential  that:  (a)  vanishes  at  a  separation  near  the  nearest 
neighbor  distance,  which  distance  depends  on  the  particular  material;  (b) , 
reproduces  the  Thcmas-Femi  potential  for  small  separations ;  (c) ,  has  for  its 
corresponding  differential  cross  section  a  simple,  easy-to-evaluate  fora;  (d) , 
has  sinple  forms  for  the  related  quantities  Na (E) ,  NS(E)  and  NW(E) ;  and  (e) , 
is  flexible  enough  to  roughly  match  any  of  the  phenomenologically  determined 
interatomic  potentials. 


& 
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We  could,  of  course,  smply  create  a  flexible  potential  form  that 
possessed  properties  (a) ,  (b)  and  (e) ,  but  there  is  no  guarantee  that  it  would 
possess  properties  (c)  and  (d) .  When  using  integral  equation  methods  for 
calculating  range  and  damage  characteristics  of  heavy  ions,  it  is  the  cross 
section,  and  not  the  potential,  that  is  used.  So  instead  of  choosing  a 
flexible  potential  form,  we  will  choose  a  flexible  cross  section  form  that 
possesses  properties  (c)  and  (d) .  Property  (b)  is  satisfied  by  requiring  that 
our  new  cross  sections  look  like  the  INS  cross  section  for  large  values  of  t. 
Because  we  can  calculate  what  the  potential  is  that  corresponds  to  a  given 
cross  section,  we  can  build  sufficient  flexibility  into  our  new  cross  sections 
in  order  to  satisfy  (e) .  Finally,  because  we  want  the  total  cross  sections  to 
be  finite,  we  will  require  that  the  kernel  f  (<p)  of  our  cross  section  fall  off 
at  least  as  feist  as  <p& ,  6>1,  as  <p  goes  to  0;  the  corresponding  potentials  will 
be  of  finite  range. 

We  express  out  new  differential  cross  sections  in  the  form 


da  »  -  ir  a2  <p~2  g  (<p)  dfi  , 


where  g  (0)  has  the  form 


g(0)  = 


<3uW>) /  0>0* 


We  will  choose  a  form  for  g u(0)  that  reproduces  the  Thamas-Fermi  cross 
section,  within  the  approximations  of  the  Lindhard  method.  For  gn,m(0) 
will  use  a  number  of  adjustable  forms.  At  this  point,  the  transition  value  0* 
is  another  free  parameter.  We  take 

gu(0)  =*  (b+0/2)/(c+0(d+0) ) ,  b=0.0074,  c=0.0376,  d=0.83.  (24 

In  Fig.  2  we  shew  a  cotparison  of  this  gu(0)  with  the  LNS  fL(0) ;  the  values  of 
ft/0)  are  taken  from  the  more  complete  tables  of  Manning,22  rather  than  from 
INS.  We  see  that  the  agreement  is  quite  good  for  0>O . 006 .  In  Fig.  3  we  show 
the  percentage  error  in  gy(0)  and  the  Winterbon,  Sigmund  and  Sanders  (WSS)  fit 
(Eq.  (21))  to  the  INS  f(0).  On  the  whole,  gu(0)  fits  fL(0)  as  well  as  does 
fy(<p)  over  the  range  0>O.OO6.  Differences  of  4-6%  are  of  no  consequence 


because  interatomic  potentials  are  not  that  well  known.  In  addition,  the  LNS 
approximations,  with  respect  to  the  exact  results,  introduce  errors  of  more 


than  4-6%. 


Par  gh,m(^)  we  use  the  form 


9n,m(<t>)  »  gn,m(*'**)  -  A  <^4  [  B  +  ]“*»  . 


Die  values  of  A  and  B  are  fixed  by  the  requirement  that  the  values  and 
derivatives  of  and  are  equal  at  0=0*.  Specifically,  we  have 


B  -  0*»V4  g.(0*)  t  (ran/40*)  ^(0*)  -  g'(0*)]-1  , 


A  -  4  (  B  +  0*^4  )n  gu(0*)  , 


0.0  L- 
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Fig.  2  —  The  kernel  of  the  differential  cross  sections  of  LNS, 
WSS,  and  the  present  work. 
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FITTING  EKXZDURE  AND  EXAMPLES 

In  Appendix  B  we  outline  a  procedure  for  obtaining  the  potential  or 
screening  function  corresponding  to  a  given  differential  cross  section,  within 
the  spirit  of  the  IMS  approximations .  By  using  this  procedure,  we  plot,  in 
Figs.  5-7,  the  screening  functions  associated  with  our  three  families  of  cross 
sections.  For  a  given  (n,m)  pair  (see  Eq.  (25) ) ,  we  show  a  family  of 
screening  factors  with  each  member  labeled  by  the  value  of  <p*. 

Given  a  phenomenological  potential,  our  procedure  for  obtaining  a 
differential  cross  section  is  as  follows:  Plot  the  screening  factor  for  the 
phenomenological  potential,  using  the  Firsov  choice  of  screening  radius, 


!<e/J)X]  NOUDNfld  DNIN33U3S 


unless  seme  specific  screening  radius  is  associated  with  the  potential. 

Ocnpare  this  plot  with  Figs.  5-7  to  obtain  the  family  member  that  best  matches 
the  given  screening  factor,  recalling  in  the  process  that  the  fit  need  not  be 
good  beyond  distances  at  which  the  potential  is  too  weak  (1-10  eV)  to  affect 
radiation  damage  and  range  calculations.  Having  obtained  a  fit  to  the 
screening  factor  with  parameters  (n,m,  0*) ,  Eq.  (23-27)  automatically  yield 
the  corresponding  cross  section. 

As  an  example,  we  consider  the  Genthcn  potential  for  copper.24  This 
potential  is  of  Thcmas-Fermi  form  out  to  0.51  A,  with  the  Lindhard  screening 
radius,  and  of  Bom-Mayer  (experiential)  form  beyond  that.  Genthon  set  the 


Fig.  6  —  The  screening  functions  corresponding  to  the  (2,3)  family 

of  cross  sections. 
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various  potentials  and  their  cross  sections,  we  plot,  in  Fig.  9,  the  Genthon 
Cu-Cu  potential,  the  potential  associated  with  our  (2,3,  .02)  cross  section, 
two  forms  of  the  Moliere  potential,  and  the  Thcmas-Fermi  potential.  The 
Moli&re  interaction  in  pictured  both  with  the  INS  screening  radius  (0.1078  A) 
that  Genthon  used,  and  with  the  Torrens  and  Robinson  Value  (0.0738  A)  ,16 
(Torrens  and  Robinson  chose  the  value  of  the  screening  radius  for  which  the 
Moli&re  potential  has  the  same  magnitude  at  the  nearest  neighbor  distance  as 
does  a  BorrHMayer  potential  whose  parameters  were  determined  from  elastic 
aonstant  data) . 


INTERATOMIC  POTENTIAL 
(V(r)  (eV),  PSEUDOLOGARITHMIC  SCALE] 


INTERATOMIC  SEPARATION  [r  (A)] 


Fig.  9  —  Comparison  of  the  Genthon  Cu-Cu  potential  with  the  Thomas- 
Fermi,  the  (2,3,  .02),  and  the  Moliere  potentials.  The  Moliere^ 
potential  is  shewn  using  both  the  ms  screening  lengthy  (0.1078  A) 
and  the  Torrens  and  Robinson  screening  length  (0.0738  A) . 


"•"Ci 


In  Fig.  10  we  shew  the  stopping  power  NS(E)  (see  Eqs.  (A2,  A5,  A13,  A15) ) 
far  copper  in  copper  based  an  each  of  the  potentials  shewn  in  the  previous 
figure.  Hie  values  for  the  Genthon  potential  were  obtained  by  the  same 
procedure23*19  as  used  by  Mueller  to  obtain  the  Moliere  results.  We  see  that 
the  properties  of  the  cross  sections  reflect  the  properties  of  the 
interactions.  For  higher  energies,  all  of  the  cross  sections  give  the  same 
results  with  the  exception  of  the  Moliere  with  the  Torrens  and  Robinson 
screening  radius,  which  is  considerably  lower.  For  lower  energies  the 
Tarrens-Rctoinscn-Mol iere  and  Genthon  results  agree,  because  both  potentials 


were  fit  to  (lew  energy)  crystalline  properties  of  copper.  The  Thamas-Fermi 
results  are  much  too  large  at  lew  energies.  Our  cross  section,  the 
(2,3,  .02) ,  gives  an  adequate  approximation  to  the  Genthcn  over  the  whole 
range  of  energies.  The  results  of  using  the  Moliere  cross  section  with  the 
Genthcn  (IDS)  screening  length  are  only  fair. 

As  further  examples,  in  Fig.  11  we  present  the  screening  factors,  and  our 
fits,  far  a  number  of  interatomic  potentials:  the  uranium-uranium  potential  of 
Genthcn;24  the  helium-tantalum  potential  of  Wilson  and  Johnson;26  and  the 
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carbon-carbon  potential  of  Genthon.  Also  pictured  are  the  Thcmas-Fermi  and 
Moliere  screening  factors.  While  our  fits  can  be  somewhat  rough,  they 
represent  their  respective  phenomenological  potentials  considerably  better 
than  does  the  Thcrnas-Fermi  screening  factor.  Consequently,  our  cross  sections 
are  an  improvement  over  the  INS  or  WSS  (Thcrnas-Fermi)  cross  sections. 

In  Figs.  12  and  13  we  show  the  stopping  power  for  C-C  and  U-U, 
respectively.  For  carbon,  as  we  might  expect  from  Fig.  11,  the  Moliere  cross 
section  overestimates  the  stopping  power  nearly  as  ouch  as  does  the  Thcrnas- 
Fermi,  but  because  the  characteristic  energy  (Z2e2/a)  is  lew  for  carbon,  even 
the  Thcrnas-Fermi  cross  section  is  reasonably  accurate.  For  uranium,  both  our 
cross  section  and  the  Moliere  give  good  results.  The  Thcrnas-Fermi  gives  much 
too  high  a  stopping  power,  and  does  so  over  a  considerable  energy  range. 


ENERGY  [E  (eV)] 


Fig.  12  —  Comparison  of  stepping  powers  for  carbon  in  carbon  based 
on  the  Genthon,  Moliere,  Thcrnas-Fermi,  and  (4,4,  .06)  cross 
sections. 


As  a  fined,  example,  we  note  that  the  phenomenological  interatomic 
potentials  are  sanetimes  of  pure  Bom-Mayer  form  V(r)  =  V'  exp  (-0'r) ,  with 
the  parameters  adjusted  so  that  the  potential  reproduces  sane  crystalline 
properties  of  the  material.27  In  these  cases  the  potential  is  clearly  not 
intended  for  use  at  high  energies;  Torrens  suggests  a  restriction  to  the  range 
far  which  V(r)<W  •  The  results  of  using  our  procedure  with  a  Bom-Mayer 
potential  is  that  we  automatically  extend  its  useful  range  for 
radiation  damage  calculations.  As  an  example,  in  Fig.  14  we  shew  the 
Gibson  #2  potential  for  copper.1  Also  shown  is  the  potential  corresponding  to 
our  (2,3,  .015)  cross  section.  Our  procedure  has  not  only  yielded  a  cross 


section  corresponding  to  the  given  potential,  but  a  cross  section  that  has  a 
wider  applicability  than  the  original  potential  could  justify. 

!v! 


INTERATOMIC  POTENTIAL  [V(r)  (eV)] 


Fig.  14  —  Ccnparison  of  the  potential  based  on  the  (2,3,  .015) 
cross  section  to  the  copper-copper  potential  of  Gibson,  et  al. 
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One  question  remains.  What  cross  section  do  we  use  when  we  have  no 
phenanenological  potential  to  guide  us?  In  the  light  of  the  wide  use  of  the 
Moliere  potential  for  radiation  damage  problems,  and  the  fact  that  the  Moliere 
fares  reasonably  well  in  most  of  the  ccuparisons  we  have  made  in  this  paper, 
we  suggest  cur  (finite  range)  version  of  the  Moliere  cross  section  with  the 
Firsov  screening  radius  as  a  universal  cross  section.  Following  the 
procedure  discussed  earlier,  choose  the  (3,2,  .019)  cross  section  to 
represent  the  Moliere  cross  section.  In  Fig.  15  we  compare  the  reduced 
stopping  power  (NS(E)/(N  it  a2  E^)  for  the  Moliere  and  (3,2,  .019)  cross 
sections.  They  are  in  good  agreement  for  all  but  the  lowest  values  of 
(reduced)  energy  transfer.  No  universal  cross  section  is  going  to  be  more 
than  adequate,  in  general,  but  in  every  case  we  have  examined  in  detail,  the 
(3,2,  .019)  cross  section  would  be  an  inprovement  over  the  Thcmas-Fermi  (INS) 


DISCUSSION 

Me  had  as  cur  goal  in  this  paper  to  create  a  class  of  cross  sections  for 
atonic  scattering  that  Mere  both  phenomenologically  based  and  easy  to  use.  We 
have  compared  our  cross  sections  to  others  in  two  ways;  by  ccrparing  the 
corresponding  potentials  or  screening  factors,  and  by  comparing  the  stepping 
powers  derived  from  each  cross  section.  We  have  not  ccnpared  total  cross 
sections  because  they  are  infinite  for  the  Thcmas-Fermi  and  Moliere  forms. 
Also,  we  have  not  ccnpared  the  square  fluctuations  in  energy  loss  (Eq.  (A4) ) 
because  tlus  quantity  depends  on  the  higher  energy  portion  of  the  kernels  of 
the  cross  sec  cions,  and  at  higher  energies  the  kernels  all  agree. 

We  can  new  ask  what  differences  might  result  from  the  use  of  these  new 
cross  sections,  relative  to  using  those  based  strictly  on  the  Thcmas-Fermi 
potential.  We  have  no  definite  answer,  but  we  can  learn  something  of  this 
matter  by  examining  the  curves  fer  stcppi.  t  power  that  we  showed  earlier.  In 
particular,  consider  Fig.  10  for  copper.  The  calculated  range  of  a  copper  ion 
with  an  initial  energy  of  more  than  a  few  keV  will  be  the  same  whether  it  is 
estimated  using  our  new  cross  sections  or  the  Thcmas-Fermi,  because  the 
stepping  powers  for  the  two  cases  are  the  same.  Similarly,  the  (initial) 
energy  deposition  curve  (energy  lest  versus  depth)  for  the  impinging  ion  will 
be  the  same  for  both  cross  sections. 

Differences  may  arise,  however,  when  the  motion  of  the  primary  knock-on 
atoms  (PKAs)  is  considered.  A  typical  PKA  will  have  an  energy  of  a  few  keV, 
and  secondary  knock-on  atoms  will  typically  have  energies  of  a  w  hundred 
electron  volts.  We  see  from  Fig.  10  that  at  these  energies  the  Thcmas-Fermi 
stepping  power  is  considerably  higher  than  is  the  phenomenologically  based 
value.  Consequently,  the  knock-on  atoms  will  have  longer  ranges  than  the 
Thcmas-Fermi  results  would  indicate,  which  may  affect  the  final  energy 
deposition  and  displacement  damage  results. 

Another  difference  arises  in  calculations  of  energy  partitioning.  An 
impinging  ion's  energy  in  given  up  to  two  mechanisms:  transfers  of  energy  to 
electrons  in  the  bulk  material,  called  inelastic  or  electronic  losses;  and 
transfers  of  energy  to  the  lattice  ions,  called  elastic  losses  or  damage 
energy.  Some  of  the  elastic  losses  produce  displacements.  If  the  elastic 
scattering  cross  section  is  lowered,  then  a  higher  proportion  of  the  initial 
energy  goes  into  inelastic  losses  and  fewer  displacements  are  produced.28'29 


APPENDIX  A:  TOTAL  CROSS  SECTIONS  AND  RELATED  QUANTITIES 
Fran  Eqs.  (9,  17-20)  we  can  write 

N  <r(E)  =  N*a2  1(6) 

N  S(E)  -  N)ra2  r  ^  e"1  J(6) 

N  a(E)  -  Mra2  [r  EL  e"1]2  K(e) 


(A1 

(A2 

(A3 


to  represent  the  total  cross  section,  stepping  power,  and  square  fluctuation 
in  energy  loss,  where 
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For  the  function  1(e)  we  can  write 

I<(0,  e  <  <p* 

!>(«)»  €><*>*, 


(A4 


(A5 


(A6 


(A7 


and  similarly  for  J(e)  and  K(e).  Recall  that  f(0)  is  of  the  form  (23-25)  and 
that  the  parameters  A  and  B  are  given  by  Eqs.  (26-27) . 

For  our  three  families  of  cross  sections  and  for  e«p* ,  we  have  the 
following  expressions: 


(n,m)  =  (3,2)  ,  y  =  B_1  e1/2  ,  (AS 

I  (e)  =  A  B"2  (2y+y2)  d+y)"2  , 


(A9 


(A10 


J<(«)  -  A  B2  [12  ln(l+y)  +  (1+y)'2  (y4-4y3-18y2-12y)]  , 

K,(€)  =  (AB6/15)  [  840  In  (1+y)  + 

(1+y)-2  (5y8-8y7+l4y6-28y5+70y4-280y3-1260y2-840y)  ]  ; 


(n,m)  =  (2,3)  ,  y  =  B_1/3  c3^4  , 


x  l€)  =  r  — xL. _ 1  1„  Il+Yli  +  -1-  tan-1  i3_Y  -| 

bt/j  3(i+y  )  •“*  l-y+y  7  2  y 


j.w  -  4»bV3  [ 

12 (1+y  )  l-y+y  7  2-y 

^(6)  =*  g  AB4  [60  In  (1+y3) 

+  (1+y3) -1  (-60y3  -30y6  +10y9  -5yU  +3y15)  ]  ; 


(n,m)  -  (4,4)  ,  y  ■  e/B  , 

X<(€)  =  1/3  A  B_1  y3  (1+y) *3  , 

J,(€)  =  ^  AB  [  (1+y)"3  (12y  +30y2  +22y3  +3y4)  -  12  ln(l+y)  ] 
K<(e)  =  A  B3  y2  (  \  y4  -y3)  (1+y)"3  +  5  b2  J.(e)  . 
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Fear  €>0  ,  for  all  of  the  families,  we  have 


I^e)  -  I <(<p*)  +(4c2)-1  [(c-2bd)(2  lny  -  P^y))  -  4bc/y 

2 
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J <(0*)  +  4  C  (4b-d)  P2(y)  +  P^y)  ]  *  , 

y=*P 


^(e)  -  K,  (<*>*)  +  J  [  y2  +(4b-2d)y  +(d2-c-2bd)  P-^y) 


+(2b(d2-2c)  -d(d2-3c) )  P2(y)  ]  * 


where 


P1(y)  -  ln[cty(d+y)  ] 

P2(y)  *  Pq”1  ln[(d+-2y-PQ)/(d+2y+-P0)]  , 


PQ  *  (d2  -  4c)1/2  . 
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APPENDIX  B:  INVERSION  MEIHOD 


Given  the  approximations  that  allowed  us  to  take  Eq.  (5)  into  Eq.  (16) , 
cne  can  easily18  obtain  the  potential  corresponding  to  a  given  cross  section 
according  to  the  prescription 


Z,  Z„  e 


dp  <P( P)  CP2  -  r2]_1/2  . 


Just  as  Eq.  (14)  is  the  first  in  a  series  of  approximations  to  the  scattering 
angle,  Eq.  (Bl)  is  the  first  tens  in  what  is  known  as  the  iirpact  expansion 
method  of  expressing  the  potential  in  terms  of  the  reduced  scattering 
angle. 81/ 

Let  us  suppose  that  we  have  a  potential  V(r)  that  vanishes  for  r>R,  and 
let  us  write  the  impact  parameter  and  R  in  reduced  the  notation 


q-P/a,  Q - R / a  . 


Then  Eqs.  (16)  and  (Bl) ,  with  the  use  of  the  approximation  (15) ,  can  be 
written 


0(q)  =  ~  2  q  dx  [X2  -  q2f 1/2  jjL  X-1  C(x) 

Jq 


dq  <P( q)  [q2  -  x2]" 1/2 


Within  the  framework  of  this  inversion  method,  the  relationship  between 
the  reduced  impact  parameter  and  the  reduced  energy  transfer  is 


q2  -  ds  s’2  g(s)  ,  0  =[TE/(rE.2)] 


w<\ 


mv/f  wwvvr 


The  cutoff  radius  is  related  to  the  differential  cross  section  by 


Q2  *  a*  <P~2  g  (0)  . 
.  o 


A  change  of  variable  frcm  q  to  <p  converts  Bq.  (B4)  to  the  form 


coo  -*  ;>  to2- *2f^2 


where  <p  and  q  are  related  by  Eq.  (B5)  and  0*  is  that  value  of  0  that  satisfies 


d 0  0~2  g(0) 


As  a  practical  matter,  we  choose  a  set  of  values  of  <p^  and  then  use  Eqs.  (B8) 
and  (B7)  to  determine  the  corresponding  values  of  the  argument  x  and  the 
screening  factor  C(x) .  in  effect,  we  form  a  table  of  C(x)  versus  x. 

As  a  example  of  the  use  of  this  procedure,  we  use  Eq.  (B7)  to  find  the 
potential  corresponding  to  fw,  the  WSS  form  of  the  INS  kernel.  In  Fig.  16  we 
shew  the  actual  Thanas-Fermi13  screening  factor  and  that  obtained  by  the  above 
procedure  using  %(0)  •  The  difference  at  large  separations  reflects  the  fact 
that  fw(0)  falls  off  more  slcwly  as  0  goes  to  0  than  would  the  equivalent 
Thomas-Fermi  kernel;  one  can  see  this  be  comparing  the  WSS  and  INS  kernels  in 
Fig.  4. 
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